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Abstract 

Considerable  effort  has  been  devoted  to  the  modeling  of  proton  exchange  membrane  fuel  cells  (PEMFCs)  as  well  as  fuel  processing  units  (FPUs). 
Many  of  these  models  consider  only  steady  state  analysis;  the  available  dynamic  models  typically  operate  only  in  simple  open  loop  configurations. 
However,  a  liquid  fuel  processor/PEMFC  stack  power  unit  for  vehicular  application  will  require  tight  integration  and  regulation  of  multiple  units 
in  order  to  function  economically  and  reliably.  Moreover,  vehicular  operation  is  inherently  dynamic  in  nature,  so  traditional  steady  state  process 
design  approaches  will  be  of  limited  value. 

This  work  addresses  a  minimum  set  of  subcomponents  necessary  for  modeling  an  overall  vehicular  power  system.  Additionally,  the  integration 
and  control  of  these  sub-units  is  addressed  so  that  the  unit  can  be  operated  as  needed  in  a  vehicular  application  by  following  a  reference  power 
trajectory.  A  number  of  design  and  operational  parameters  can  be  adjusted  and  the  impact  on  system  performance  studied.  Based  on  this  preliminary 
analysis,  heuristics  are  developed  for  optimal  operation  and  design. 

©  2006  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

Dwindling  supplies  of  fossil  fuels  have  caused  researchers 
and  policy  makers  to  explore  alternative  energy  sources.  As  of 
2003,  energy  consumption  by  the  transportation  sector  com¬ 
prised  27%  of  total  US  energy  consumption  and  67%  of  all 
US  petroleum  consumed  [1].  Consequently,  one  of  the  largest 
alternative  energy  research  initiatives  lies  in  transportation  tech¬ 
nology.  Developing  technologies  that  reduce  the  demand  on 
petroleum  products  has  environmental  benefits  from  reduced 
consumption  in  addition  to  economic  and  political  benefits 
from  the  reduced  dependence  on  foreign  oil  imports.  Despite 
continued  improvements  in  engine  and  fuel  injection  technol¬ 
ogy,  which  have  increased  the  average  fuel  economy  of  pas¬ 
senger  vehicles  from  13.5  miles  per  gallon  (mpg)  in  1970  to 
22.1  mpg  in  2002,  both  the  total  number  of  vehicles  and  the 
average  annual  distance  travelled  per  vehicle  have  steadily 
increased  over  the  same  time  period.  The  net  result  is  an 
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approximately  70%  increase  in  the  amount  of  motor  fuel 
consumed  [1]. 

Even  considering  improved  fuel  economy,  current  internal 
combustion  engines  (ICEs)  achieve  an  energy  efficiency  of 
approximately  20%  [48].  Further  efficiency  increases  are  un¬ 
likely  without  significant  modifications  and  technological  break¬ 
throughs.  Therefore,  a  significant  amount  of  research  has  been 
directed  towards  alternative  fuel  technologies,  such  as  hybrids, 
plug-in  hybrids,  and  fuel  cell  vehicles  (FCVs).  Hybrids  combine 
a  small  ICE  with  a  battery  system  connected  to  an  electric  mo¬ 
tor/generator.  While  the  vehicle  is  powered  largely  by  the  ICE, 
the  electric  motor  provides  supplemental  power  from  the  battery 
pack  during  high-load  conditions  such  as  acceleration.  During 
lower  load  conditions  such  as  idling  or  cruising,  surplus  power 
produced  by  the  ICE  drives  the  electric  motor  which  is  then  oper¬ 
ated  in  reverse  as  a  generator  to  recharge  the  batteries.  Likewise, 
regenerative  braking  techniques  can  be  employed  during  decel¬ 
eration  to  recover  some  of  the  kinetic  energy  of  the  vehicle, 
converting  it  back  to  electrical  energy  stored  in  the  batteries. 
Since  these  vehicles  incorporate  already  existing  ICE  technol¬ 
ogy,  they  have  reached  a  more  advanced  stage  of  development 
with  a  number  of  models  such  as  the  Toyota  Prius,  Honda  In- 
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Nomenclature 

Amem  PEMFC  stack  membrane  area  for  each  cell  (cm2) 
c  total  molar  concentration  of  gaseous  species  in 
methanol  reformer  (mol  m-3) 

Ci  molar  concentration  of  species  i  in  methanol  re¬ 
former  (mol  m-3) 

Cq2  saturation  concentration  of  oxygen  in  the  PEMFC 

membrane  (mol  m3) 

C  capacitance  of  PEMFC  stack  (F) 

CaUx  auxiliary  power  capacity  (A  h  at  120  V) 

Cp  constant  pressure  heat  capacity  (Jmol-1  or 
Jkg-1) 

Cp?Cat  Heat  capacity  of  catalyst  bed  in  reformer  (J  kg- 1 ) 
Fp,gas  Heat  capacity  of  reformer  gas  (J  mol-1) 

Es  energy  in  PEMFC  stack  body  (J) 

Es  rate  of  change  of  energy  in  PEMFC  stack  body 
(J  min-1) 

Eq j  Equilibrium  driving  force  of  reaction  j  in 
methanol  reformer 
T  Faraday  constant  (C  mol-1) 
h  channel  height  in  cross-flow  heat  exchanger  (m) 
(hA)  lumped  convective  heat  transfer  rate  (J  min-1  K) 
Hi  specific  enthalpy  of  species  i  (Jkg-1) 

A  Hrxn  ?  j  heat  of  reaction  j  in  methanol  reformer  ( J  mol  - 1 ) 
i  current  density  in  PEMFC  stack  (A  cm-2) 

I  current  in  PEMFC  stack  (A) 
kj  rate  constant  for  reaction  j  in  methanol  reformer 
(molm-3  min-1  kPa) 

kY  valve  proportionality  constant  (m3  min-1  atm) 
Kmx  auxiliary  power  recharging  proportionality  con¬ 
stant  (W  A-1  h) 

Kc  PI  controller  gain 

Keq  j  equilibrium  constant  of  reaction  j  in  methanol  re¬ 
former  (kPax) 

Lx  length  of  cold  channel  in  cross-flow  heat  ex¬ 
changer 

Ly  length  of  hot  channel  in  cross-flow  heat  exchanger 

nti  mass  of  species  i  (kg) 

riii  mass  flow  rate  of  species  i  (kg  min-1) 

Mi  molecular  weight  of  species  i  (kg  mol-1) 

N  number  of  cells  in  PEMFC  stack 
Pan  pressure  in  anode  channel  of  PEMFC  stack  (atm) 
Pea  pressure  in  cathode  channel  of  PEMFC  stack 
(atm) 

Pdesired  desired  system  power  level  (W) 

Pi  partial  pressure  of  species  i  in  the  reformer  (kPa) 
Psp  PEMFC  stack  power  setpoint  (W) 

Pstack  PEMFC  stack  power  (W) 

A  P  pressure  drop  across  a  valve  (atm) 
q  heat  flux  (J  m-2  s) 

•  t 

Gcool  cooling  rate  on  PEMFC  stack  (J  min-1) 

r  j  reaction  rate  of  reaction  j  in  methanol  reformer 
(mol  m  3  min) 

R  ideal  gas  constant  (J  mol-1  K) 


Pact  cell  resistance  due  to  oxygen  activation  losses 
(ohm) 

Pohm  cell  resistance  due  to  transport  (Ohmic)  losses 
(ohm) 

Rr  radius  of  reformer  reactor  (m) 
s  Laplace  domain  variable 
t  time  (min) 

T  methanol  reformer  temperature  (K) 

Pwaii  external  wall  temperature  of  methanol  reformer 
(K) 

vz  gas  velocity  in  axial  direction  of  methanol  re¬ 
former  (m  min- 1 ) 

U  convective  heat  transfer  coefficient  (J  m-2  s  K) 

V  volumetric  flow  rate  (m3  min-1) 

Fact  Dynamic  activation  overvoltage  (V) 

Vceii  voltage  of  a  single  cell  within  the  PEMFC  stack 

(V) 

Fstack  PEMFC  stack  voltage  (V) 

z  axial  dimension  of  methanol  reformer  (m) 

Greek  symbols 

y  overall  volume  fraction  of  methanol  reformer 
rx  weighting  factor  for  objective  x  in  multiobjective 
optimal  design  framework 
6  porosity  of  catalyst  pellets  in  methanol  reformer 

6bed  void  fraction  of  packed  catalyst  bed  in  methanol 
reformer 

77 act  steady  state  activation  overpotential  in  PEMFC 
cell  (V) 

77eff  reaction  effectiveness  factor  in  methanol  reformer 
Vi  stoichiometric  coefficient  of  species  i  in  PEMFC 
stack  reactions 

emperical  modeling  coefficient  in  PEMFC  stack 
model 

p  density  (kg  m- 3 ) 

yOcat  bed  density  of  packed  catalyst  pellets  (kg  m-3) 
tf  First  order  filter  time  constant  (min) 
x\  PI  controller  integral  time  constant  (min) 
tpx  cost  of  objective  x  in  multiobjective  optimal  de¬ 
sign  formulation 

O  multiobjective  optimal  design  cost  function 


sight,  Honda  Civic  hybrid,  Honda  Accord  hybrid,  and  the  Ford 
Escape  hybrid  already  commercially  available  and  rapidly  gain¬ 
ing  popularity.  Fuel  economies  range  from  29mpg  (city)  for 
the  high-power  Honda  Accord  hybrid  to  66  mpg  (highway)  for 
the  Honda  Insight,  according  to  EPA  estimates.  Regardless  of 
the  model,  fuel  economy  is  significantly  higher  than  the  current 
national  average. 

Plug-in  hybrids  are  similar  in  many  respects  to  conventional 
hybrids,  except  they  feature  a  smaller  ICE  with  a  larger  electric 
motor/generator  and  battery  pack.  As  the  name  implies,  plug-in 
hybrids  will  need  to  connect  to  an  external  electrical  source  to 
charge  the  battery  pack  completely.  When  fully  charged,  plug-in 
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hybrids  will  operate  entirely  on  electrical  power  until  the  power 
level  in  the  battery  pack  drops  to  a  certain  level.  It  is  expected 
that  they  will  be  able  to  drive  20-35  miles  on  batteries  alone  be¬ 
fore  the  ICE  is  required.  Once  the  ICE  starts,  operation  will  be 
more  similar  to  a  traditional  hybrid  with  power  control  system 
only  recharging  the  batteries  enough  from  the  ICE  to  maintain 
hybrid  operation.  Once  the  operator  reaches  a  destination  with 
proper  charging  equipment,  the  vehicle  can  be  plugged-in  and 
the  batteries  fully  recharged.  Given  standard  driving  habits,  this 
operational  paradigm  could  result  in  a  large  percentage  of  trans¬ 
portation  energy  coming  from  electricity  rather  than  petroleum 
fuels,  in  addition  to  the  higher  fuel  economy  of  hybrid  elec¬ 
tric  vehicles.  Some  estimates  have  placed  the  reduction  in  oil 
for  transportation  by  as  much  as  74%  -  assuming  full  market 
penetration  -  by  switching  to  the  plug-in  hybrid  configuration 
[48]. 

Both  traditional  hybrids  and  plug-in  hybrids  offer  signifi¬ 
cant  improvements  over  traditional  ICE  vehicles  in  terms  of 
fuel  economy,  but  still  require  gasoline  or  diesel  fuel.  In  order 
to  entirely  eliminate  the  need  for  fossil  fuels,  another  prominent 
research  area  involves  the  use  of  hydrogen  as  fuel.  Although 
it  is  possible  to  combust  hydrogen  in  a  traditional  ICE,  signif¬ 
icant  effort  has  been  devoted  to  fuel  cells,  particularly  proton 
exchange  membrane  fuel  cells  (PEMFCs),  which  oxidize  hy¬ 
drogen  electrochemically  to  produce  electrical  power.  Due  to 
the  lower  operating  temperature  of  PEMFCs,  it  is  believed  that 
the  energy  efficiency  of  hydrogen  fuel  cell  powered  vehicles  can 
exceed  40%.  While  the  ideal  fuel  source  for  fuel  cells  is  pure 
hydrogen,  there  are  a  number  of  technological,  economic,  and 
infrastructure-related  hurdles  that  make  it  unlikely  that  the  first 
commercially-available  FCVs  will  run  directly  on  compressed 
gaseous  or  liquid  hydrogen.  Instead,  early  FCVs  will  likely  be  fu¬ 
eled  by  liquid  hydrocarbons  such  as  methanol  (MeOH),  ethanol 
(EtOH),  gasoline,  or  naphtha,  which  is  then  converted  on-board 
into  a  hydrogen-rich  fuel  stream  [49].  While  it  would  be  con¬ 
venient  to  continue  to  use  gasoline  as  a  hydrogen  source  for 
FCVs  [18],  environmental  concerns  over  greenhouse  gas  emis¬ 
sions  will  likely  dictate  a  switch  to  lighter  hydrocarbons  such  as 
methanol  or  ethanol  which  offer  a  higher  ratio  of  hydrogen  to 
carbon  dioxide.  In  fact,  the  reforming  of  methanol  can  yield  a 
hydrogen  to  carbon  dioxide  ratio  as  high  as  3: 1,  whereas  the  par¬ 
tial  oxidation  of  octane  (as  an  approximation  for  gasoline)  can 
only  achieve  a  maximum  ratio  of  9:8.  Due  to  the  relatively  high 
hydrogen/carbon  ratio  of  methanol,  several  groups  have  begun 
studying  methanol  reforming  specifically  for  the  production  of 
hydrogen  for  PEMFC  applications  [7,27,30]. 

Considerable  work  has  already  been  done  to  develop  FCVs 
running  directly  on  pure  hydrogen.  GM,  Daimler-Chrysler, 
BMW,  and  others  have  working  prototypes.  However,  the  re¬ 
search  and  design  is  at  a  much  earlier  stage  for  vehicles  that  fea¬ 
ture  on-board  hydrogen  production  from  liquid  fuels.  Overall, 
the  systems  will  be  considerably  more  complicated,  since  the 
equipment  needed  to  convert  the  liquid  fuel  to  hydrogen-rich 
gas  and  likely  purify  it  has  its  own  set  of  energy  demands  and 
dynamic  behaviors.  Not  surprisingly,  a  number  of  design  con¬ 
siderations  will  impact  the  cost  and  efficiency  of  such  a  plant  as 
well  as  its  operability  and  controllability.  One  major  concern  is 


start-up  time:  how  long  will  it  take  for  a  reformer  to  heat  to  oper¬ 
ating  temperature  and  begin  producing  enough  hydrogen  to  run 
the  fuel  cell?  The  cold  start  of  methanol  reformers  is  currently 
under  investigation  and  strategies  are  being  developed  to  mini¬ 
mize  warm-up  time  [28,15].  Nevertheless,  in  order  to  “turn  the 
key  and  go”,  FCVs  will  likely  require  battery  networks  similar 
to  today’s  gasoline-electric  hybrids  in  order  to  power  the  vehi¬ 
cle  while  the  hydrogen  production  and  fuel  cell  systems  come 
online. 

Additionally,  heat  integration  will  be  crucial  to  the  overall  ef¬ 
ficiency  of  the  vehicle.  The  fuel  cell  itself  generates  excess  heat, 
and  depending  on  the  type  of  liquid  fuel  conversion,  that  process 
may  be  at  elevated  temperature  as  well.  In  order  to  increase  effi¬ 
ciency,  some  of  this  heat  may  be  recovered  by  using  it  to  preheat 
the  feed  stream  to  the  reformer  unit  as  well  as  the  air  feed  to  the 
cathode.  However,  this  sort  of  integration  couples  the  dynamics 
of  various  units  and  could  lead  to  unforeseen  behaviors  under 
transient  driving  conditions.  Therefore,  it  is  important  to  closely 
study  the  design  and  integration  of  any  combined  liquid  fuel 
processor  and  fuel  cell  stack  system  not  only  from  a  steady  state 
perspective,  but  from  a  dynamic  one  as  well. 

The  following  work  is  divided  into  several  sections.  First,  the 
development  and  implementation  of  the  core  components  within 
the  overall  system  model  is  discussed.  This  is  followed  by  the 
assembly  of  the  components  into  a  complete  system  model  and 
a  description  of  the  necessary  control  loops  for  dynamic  power- 
level  operation.  Using  the  system  structure  presented,  design 
sensitivity  studies  were  performed  and  the  results  examined. 
Finally,  future  directions  for  this  work  are  considered. 

2.  Component  modeling 

The  liquid-fueled  PEMFC  power  plant  is  by  design  a  com¬ 
bination  of  subsystems  coupled  together.  In  order  to  study  the 
behavior  of  the  system  as  a  whole,  it  is  necessary  to  have  appro¬ 
priate  models  of  all  constituent  components.  At  the  very  mini¬ 
mum,  models  for  the  PEMFC  stack  and  a  fuel  processing  unit 
are  required.  Ultimately,  the  fuel  processing  unit  will  itself  be 
comprised  of  several  sub-units  including  reforming,  water-gas 
shift,  and  preferential  oxidation  reactors  or  reaction  zones  within 
a  single  reactor.  This  complicated  arrangement  is  necessary  to 
reduce  the  amount  of  CO  and  other  PEMFC  catalyst  poisons 
present  in  the  reformate  stream.  However,  in  an  effort  to  reduce 
the  complexity  of  the  overall  system  model,  a  single  model  of 
a  methanol  reformer  is  used  individually  as  the  fuel  process¬ 
ing  unit  (FPU).  Since  the  reforming  reaction  is  endothermic,  the 
dynamics  of  the  reformer  should  be  the  slowest  of  all  the  com¬ 
ponents  within  the  fuel  processor,  and  thus  should  be  a  good 
indicator  of  the  behavior  of  the  fuel  processing  unit  as  a  whole. 

A  model  for  the  PEMFC  stack  is  the  other  critical  compo¬ 
nent.  While  very  sophisticated  models  have  been  developed  for 
various  fuel  cells  and  stacks  [20,45,22],  their  computational  cost 
is  quite  large  and  generally  prohibitive  for  dynamic  studies,  par¬ 
ticularly  at  the  system  level.  The  desire  for  this  work  is  to  have 
a  relatively  simple  model  that  provides  a  reasonably  accurate 
account  of  the  overall  material  and  energy  flows  of  the  system 
as  well  as  the  voltage  behavior  of  the  stack  as  a  whole  in  an 
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efficient  manner.  A  suitable  model  has  been  adapted  from  the 
literature  to  accomplish  these  goals  [51,38]. 

Finally,  a  variety  of  ancillary  equipment  is  required  to  build 
an  actual  system.  This  includes  pumps,  burners/heaters,  blow¬ 
ers,  heat  exchangers/evaporators,  auxiliary  power  sources,  and 
power  controllers/converters.  However,  many  of  these  items  add 
needless  complexity  to  the  model  while  providing  minimal  ad¬ 
ditional  insight  into  the  process.  Nevertheless,  two  items  from 
this  list  bear  further  discussion.  The  reformer  and  PEMFC  stack 
operate  at  considerably  different  temperatures  and  so  it  becomes 
necessary  to  cool,  heat,  and  evaporate  various  process  streams 
at  different  points  to  operate  safely  and  efficiently.  Therefore, 
there  is  the  opportunity  for  heat  integration  to  use  one  stream  to 
heat  or  cool  a  second  in  the  process.  This  heat  integration  more 
tightly  couples  the  operation  of  all  units,  but  should  ultimately 
allow  for  a  higher  energy  efficiency  for  the  system.  Hence,  there 
is  a  need  for  some  basic  heat  exchanger  modeling  to  examine 
the  effects  of  different  heat  integration  schemes. 

In  addition  to  the  heat  exchangers,  the  role  of  the  auxiliary 
power  supply  is  also  important.  By  examining  the  amount  of 
auxiliary  energy  and  power  needed  to  operate  a  particular  sys¬ 
tem  configuration,  one  can  determine  what  type  and  size  of  aux¬ 
iliary  power  source  is  needed.  From  that  data,  economic,  size, 
and  weight  constraints  may  dictate  whether  a  particular  design 
is  practical.  A  more  detailed  treatment  of  each  component  is 
presented  below. 

2.7.  Methanol  reformer 
2.1.1.  Modeling  assumptions 

In  order  to  clarify  development  and  speed  computation,  a 
number  of  modeling  assumptions  and  simplifications  have  been 
made.  One  significant  assumption  is  that  the  reactor  behaves  as 
if  it  is  plug-flow.  Therefore,  the  only  spatial  gradients  in  this 
model  will  be  in  the  axial  direction  along  the  length  of  the  re¬ 
actor.  Given  the  relatively  large  flow  rate  and  the  packed-bed 
conditions  inside  the  reactor,  this  should  provide  a  reasonable 
approximation.  It  is  also  assumed  that  all  chemical  components 
remain  in  the  gas  phase,  so  that  the  system  can  be  modeled  as  a 
single  phase.  Furthermore,  the  gasses  are  assumed  to  obey  the 
ideal  gas  law  in  order  to  avoid  the  solution  of  nonlinear  equa¬ 
tions  of  state  at  each  step  of  integration.  This  is  a  reasonable 
assumption,  considering  the  operating  conditions  of  this  reactor 
with  temperatures  greater  than  200  °C  and  pressures  well  below 
10  atm.  A  constant  pressure  assumption  is  made  for  the  sake  of 
simplicity.  By  using  the  constant  pressure  assumption  together 
with  the  ideal  gas  law  and  an  assumption  of  no  spatial  gradi¬ 
ents  in  mass  flux,  the  superficial  velocity  can  be  identified  as  a 
function  of  position. 

Finally,  the  heat  transfer  in  the  reactor  must  be  considered. 
For  a  reactor  with  pressure  driven  flow,  it  is  reasonable  to  assume 
that  the  Peclet  number,  the  ratio  of  convective  heat  transfer  to 
conductive  heat  transfer:  pCpkv  L ,  is  quite  large,  so  conduction 
within  the  gas  and  between  the  gas  and  the  solid  catalyst  pellets 
is  neglected.  Additionally,  the  catalyst  bed  is  assumed  to  have 
the  same  temperature  distribution  as  the  gas  flowing  through  it. 


The  transfer  of  heat  into  the  reactor  from  the  surroundings  is 
treated  empirically  as  convective  heat  transfer  using  Newton’s 
law  of  cooling,  which  states  that  the  heat  flux  into  the  reactor 
is  proportional  to  the  difference  between  the  gas  temperature 
inside  the  reactor  and  the  specified  external  wall  temperature 
given  in  the  following  equation: 

q(z)  =  U (Twan  -  T(z))  (1) 

The  parameter  U  is  the  heat  transfer  coefficient  which  is  either 
determined  experimentally  or  through  an  experimental  correla¬ 
tion  such  as  the  Colburn  analogy  [16].  Values  will  typically  be  in 
the  range  of  25-100  J  m-2  s-1  K.  For  the  purpose  of  this  model, 
it  is  assumed  that  the  external  wall  temperature  is  uniform  and 
can  be  maintained  at  a  desired  temperature  by  a  heating  element. 
In  practice,  the  heat  will  likely  be  provided  by  the  catalytic  com¬ 
bustion  of  methanol  and/or  the  hydrogen  remaining  in  the  anode 
off  gas.  The  utilization  of  the  excess  hydrogen  in  the  off  gas  will 
be  part  of  the  heat  integration  process  and  critical  to  the  overall 
fuel  efficiency  of  the  vehicle. 

2.1.2.  Reaction  network  and  kinetics 

Given  the  current  interest  in  hydrogen  production  from  liquid 
fuels  such  as  methanol,  a  number  of  studies  have  recently  ap¬ 
peared  in  the  literature  characterizing  various  types  of  catalysts 
for  the  numerous  routes  of  hydrogen  production  from  methanol. 
These  studies  include  methanol  decomposition,  steam  re¬ 
forming,  partial  oxidation,  and  combined/autothermal/oxidative 
methanol  reforming — a  combination  of  steam  reforming  and 
partial  oxidation.  Additionally,  several  novel  reactor  configura¬ 
tions  have  been  proposed  with  improved  heat  and  mass  transfer 
properties  over  traditional  packed-bed  reactors. 

Of  all  the  potential  catalysts  studied,  Cu  and  mixtures 
of  Cu/ZnO  on  AI2O3  have  received  the  most  attention 
[34,47,40,39,2,3,12].  These  catalysts  have  been  in  use  commer¬ 
cially  for  several  decades  in  the  synthesis  of  methanol.  With  the 
addition  of  water  as  a  reactant,  they  work  very  well  to  catalyze 
the  steam  reformation  of  methanol  to  produce  H2  and  CO2  with 
a  high  conversion  of  methanol  and  a  high  selectivity  to  H2  and 
subsequent  low  selectivity  towards  CO.  This  is  particularly  im¬ 
portant  for  PEM  fuel  cell  applications,  since  CO  is  a  poison  to 
the  proton  exchange  catalysts  even  at  very  low  concentrations. 
One  of  the  reasons  for  the  low  selectivity  towards  CO  is  that  the 
catalyst  also  facilitates  the  water-gas  shift  reaction: 

CO  +  H20  ►  C02  +  H2  (2) 

which  is  driven  towards  the  right  with  a  stoichiometric  excess  of 
water.  However,  steam  reforming  is  a  relatively  endothermic  re¬ 
action  ( A  77rxn  =  49  kJ  mol- 1 )  which  requires  an  external  source 
of  heat,  a  definite  design  hurdle  for  a  portable  application  such 
as  a  fuel  cell  vehicle.  To  make  matters  worse,  these  catalysts 
exhibit  slow  reaction  rates  and  therefore  need  to  operate  at  tem¬ 
peratures  exceeding  250  °C.  This  combination  of  factors  makes 
traditional  packed-bed  reactors  using  these  catalysts  thermally 
sluggish,  which  in  turn  makes  it  difficult  to  respond  to  rapid 
changes  in  load. 

Others  such  as  Agrell  et  al.  [2]  have  examined  the  effects  of 
using  zirconia  (Zr02)  as  an  alternate  (or  additional)  support  ma- 
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terial  for  Cu/ZnO  with  some  performance  gains  during  steam  ref¬ 
ormation.  However,  these  catalysts  exhibited  a  higher  light-off 
temperature  and  elevated  CO  production  compared  to  the  AI2O3 
only  catalysts  when  used  in  an  oxidative  steam  reforming  en¬ 
vironment.  Moving  away  from  steam  reforming  entirely,  others 
have  discussed  direct  and  oxidative  decomposition  of  methanol 
on  Pd/Si02  (references  contained  in  [23])  and  the  use  of  in¬ 
dium  tin  oxide/alumina  (ITO/AI2O3)  for  the  partial  oxidation  of 
methanol  (POM)  [29],  shown  in  the  following  equation: 

MeOH  +  ±02  — >  2H2  +  C02  (3) 

Direct  decomposition  of  methanol  is  impractical  for  portable 
fuel  cell  purposes,  in  that  it  is  both  highly  endothermic  (A  Hrxn  = 
+90.6  kJ  mol-1)  and  it  produces  a  large  amount  of  CO  as 
a  product.  Conversely,  POM  is  highly  exothermic  (A Hrxn  = 
—  1 92  kJ  mol- 1 )  which  makes  it  difficult  to  control.  Additionally, 
the  maximum  yield  of  H2  is  only  2  mol  mol- 1  MeOH  with  POM 
compared  to  a  possible  3  mol  mol-1  MeOH  with  steam  reform¬ 
ing,  with  a  fairly  significant  production  of  CO  as  a  byproduct. 
Thus,  the  most  promising  technique  for  hydrogen  generation  for 
portable  PEM  fuel  cell  applications  is  oxidative  steam  reforming 
of  methanol  (OSRM),  also  known  as  combined  steam  reforming 
(CSR)  or  autothermal  reforming,  when  the  system  is  operated 
adiabatically.  In  addition  to  the  water  and  methanol  needed  for 
steam  reforming,  a  sub- stoichiometric  amount  of  O2  (or  air)  is 
added  to  the  feed  so  that  some  methanol  is  consumed  by  the 
POM  pathway.  The  overall  OSRM  reaction  is  given  by: 

MeOH  +  (1  -  x)H20  +  \x02  — >  (3  -  jc)H2  +  C02  (4) 

where  v  is  the  fraction  of  methanol  consumed  by  the  POM  mech¬ 
anism  [30] .  Consequently,  the  overall  heat  of  reaction  for  OSRM 
depends  on  v  and  is  given  by  A  Hrxn  =  49.5  —  241 .5x  kJ  mol-1 . 
Adjusting  reaction  conditions  so  that  v  ^  0.2  results  in  an  adia¬ 
batic  or  autothermal  system  in  which  the  overall  energy  change 
is  neutral.  There  are  numerous  advantages  to  OSRM  (see  Turco 
et  al.  [47]  and  references  therein),  (i)  The  heat  released  by  the 
POM  is  used  to  drive  the  reforming  reaction.  When  the  O2  (air) 
stream  is  regulated  properly,  no  additional  heat  will  be  need  to 
be  added  to  the  reactor  to  sustain  operation;  (ii)  By  using  the 
excess  heat  of  POM  for  steam  reforming,  the  reactor  tempera¬ 
ture  is  more  easily  controlled  than  an  exothermic  POM  reactor 
and  responds  to  transients  more  quickly  than  an  endothermic 
reformer;  (iii)  The  presence  of  water  in  the  system  drives  the 
water-gas  shift  reaction  towards  H2  and  CO2  which  helps  elim¬ 
inate  the  CO  produced  from  incomplete  oxidation  in  the  POM 
pathway;  (iv)  OSRM  has  a  slightly  lower  H2  yield  than  tradi¬ 
tional  steam  reforming,  but  the  elimination  of  the  need  for  an 
exogenous  heat  supply  and  improved  dynamic  response  makes 
up  for  this  loss.  Moreover,  the  yield  is  still  noticeably  better  than 
POM  alone. 

Unfortunately,  very  little  information  is  available  in  the  liter¬ 
ature  about  the  overall  kinetics  of  OSRM.  In  order  to  prototype 
the  design  methodology,  it  was  decided  to  focus  on  steam  re¬ 
forming,  for  which  more  quantitative  information  exists.  For 
the  present  version  of  the  model,  the  kinetics  of  the  reformation 
of  methanol  were  obtained  from  Mizsey  et  al.  [34].  The  reaction 


network  was  given  as  follows: 

MeOH  4 — »  CO  +  2H2  (+90.6  kJ  mol”1) 

methanol  decomposition  (D)  (5) 

CO  +  H20  4 — >  C02  +  H2  (-41.1  kJ  mol”1) 

water-gas  shift  (W)  (6) 


MeOH  +  H20  C02  +  3H2  (+49.5  kJ  mol”1) 

steam  reforming  (R)  (7) 


Each  of  these  reactions  has  an  experimentally-determined, 
temperature-dependent  equilibrium  constant,  given  by  the  fol¬ 
lowing  equations: 


d  =  1.718  x  1014  exp  | 
^feq,w  =  9.543  x  10-3  exp 
^eq,R  =  1.849  x  1010  exp  | 


95418\ 
RT  ) 

39876  \ 
RT  ) 

56087  \ 
RT  ) 


(8) 

(9) 

(10) 


Using  the  equilibrium  coefficients,  equilibrium  driving  force  ex¬ 
pressions  can  be  developed: 

>2 
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Finally,  the  overall  rate  expressions  for  the  three  reactions  are 
given  by: 


HD  =  ^D^MeOHEqo  (14) 

rw  =  kwPcoEqw  (15) 

rR  =  ^RPMeOHEqR  (16) 

where  the  rate  constants  kx  are  assumed  to  obey  the  Arrhe¬ 
nius  temperature  dependence  with  activation  energy  and  pre¬ 
exponential  factors  given  in  Table  1 .  It  should  be  noted  that  this 
reaction  mechanism  was  determined  with  water  in  excess.  Ad¬ 
ditionally,  this  mechanism  does  not  incorporate  the  methanol 
oxidation  pathway,  which  is  required  for  autothermal  methanol 
reforming. 


Table  1 


Arrhenius  parameters  for  methanol  decomposition  (obtained  from  Mizsey  et  al. 
[34]) 


Reaction 

Activation  energy 

Pre-exponential  factor 

(kJ  mol  - 1 ) 

(mol  g  cat- 1  s- 1  kPax) 

Methanol  decomposition 

76(±4%) 

1.12(±12%) 

Water-gas  shift 

50(±25%) 

0.0023(±32%) 

Steam  reforming 

81(±7%) 

6.75(±34%) 
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Note  that  by  adding  the  stoichiometric  coefficients  of  the 
decomposition  and  water-gas  shift  reactions,  one  obtains  the 
stoichiometry  of  the  steam  reforming  reaction.  Thus,  this  system 
is  not  actually  three  independent  reactions.  Since  only  two  of  the 
three  reactions  are  linearly  independent,  only  two  of  the  species 
reaction  rates  are  distinct  and  the  remainder  are  expressed  as  a 
combination  of  the  first  two  as  demonstrated  in  the  following 


equations: 

rco2  =  rw  +  rR  (17) 

rco  =  fD  —  ^w  (18) 

m2o  =  ~rco2  (19) 

^MeOH  =  -(^C02  +  ^CO)  (20) 

m2  =  3rco2  +  2rco  (21) 

2.1.3.  Model  equations 


The  model  equations  for  this  system  consist  of  the  stan¬ 
dard  multi-component  mass  and  energy  balances  for  cylindri¬ 
cal,  plug-flow  geometry  with  corresponding  modifications  for 
the  packed  nature  of  the  bed  and  the  form  of  heat  transfer  into 
the  reactor.  A  void  fraction  <rbed  is  specified,  which  denotes  what 
amount  of  the  bed  is  unoccupied  by  pellets.  Values  of  the  void 
fraction  vary  widely  based  on  the  size  and  geometry  of  the  pack¬ 
ing  material.  However,  for  packed-bed  catalytic  reactors,  values 
of  0. 3-0.5  are  common  [21].  The  methanol  reforming  reactor 
used  by  Dams  et  al.  [17]  is  estimated  to  be  0.37.  Additionally, 
the  ratio  of  pore  volume  to  total  pellet  volume  is  the  porosity, 
denoted  €.  By  combining  the  void  fraction  and  pellet  porosity, 
the  volume  fraction  y  available  to  the  gas  within  the  reactor  is 
given  by  the  following  equation: 

x  =  e bed  +  (l -<fbed)e  (22) 

At  this  point,  one  may  start  from  very  general  mass  and  energy 
balances  and  begin  simplifying  as  demonstrated  in  Nauman,  pp. 
531-533  [35].  However,  this  particular  analysis  has  been  per¬ 
formed  previously  for  the  methanol  reforming  system  by  Dams 
et  al.  [17],  giving  rise  to  the  component  mass  balance: 

dc[  d 

y—  +  <fbed  —  (Vz  Cj)  =  r)eif  pcat  r,  (23) 

at  oz 

and  the  energy  balance: 

dT 

’gaS~dz 

=  2— -(Twan  —  T)  —  T]e ff  pcat  ^  ^  A Hrxnj  r j  (24) 

r  je{  D,W,R} 

In  Eqs.  (23)  and  (24)  above,  c*  is  molar  concentration  of 
species  i,  while  c  is  the  molar  concentration  of  all  gaseous 
species.  Additionally,  ?7eff  is  the  effectiveness  factor  which  can  be 
used  to  account  for  mass  transfer  limitation  of  the  rate  of  reaction 
within  the  pellet;  C  p  gas  is  the  molar  heat  capacity  (J  mol- 1  K- 1 ) 
of  the  gas  phase;  Cp5Cat  is  the  specific  heat  capacity  (J  kg-1  K-1) 
of  the  catalyst  pellets;  pc at  is  the  packed  density  of  the  catalyst 
pellets;  R  is  the  radius  of  the  reactor  tube.  Finally,  A Hrxnj  is 


dT 

(ycC p,gas  +  Pc atCp,cat)“  ^bcdOz^Cp 

ot 


the  heat  of  reaction  for  the  jth  reaction  pathway  in  the  system. 
Combined  with  the  ideal  gas  law: 


RT 


(25) 


this  completes  the  modeling  of  the  packed-bed  methanol  re¬ 
former. 


2.1.4.  Implementation 

The  governing  mass  and  energy  balances  for  the  reformer 
are  partial  differential  equations  (PDEs)  with  a  single  spatial 
dimension.  Very  few  nonlinear  PDEs  have  analytical  solutions, 
so  it  becomes  necessary  to  rely  on  numerical  techniques.  A  va¬ 
riety  of  methods  exists  for  the  numerical  solutions  of  differ¬ 
ent  classes  of  PDEs.  The  reformer  model  consists  of  six  cou¬ 
pled  parabolic  PDEs  with  a  simple  linear  geometry,  so  tech¬ 
niques  such  as  collocation,  finite  differences,  and  the  method 
of  lines  would  be  suitable  solution  strategies.  The  method  of 
lines  uses  finite  differences  to  discretize  the  spatial  deriva¬ 
tives  across  a  domain  of  interest,  thereby  converting  a  system 
of  PDEs  into  a  larger  system  of  coupled  ordinary  differential 
equations  (ODEs)  depending  only  on  time  [42].  Many  robust 
algorithms  and  packages  exist  for  solving  systems  of  ODEs 
[32,33,26,41],  which  generally  makes  this  method  more  reli¬ 
able  than  a  pure  finite  difference  approach  where  both  spatial 
and  temporal  derivatives  are  discretized.  Since  the  time  deriva¬ 
tives  of  Eqs.  (23)  and  (24)  can  be  isolated  algebraically,  this  sys¬ 
tem  is  an  excellent  candidate  for  solution  through  the  method  of 
lines. 

It  was  determined  that  at  least  50  discretization  intervals  were 
required  to  sufficiently  reduce  spatial  discretization  error,  which 
results  in  a  system  of  300  coupled  ODEs  (6  states  at  50  points). 
Given  the  large  number  of  equations,  the  system  was  imple¬ 
mented  in  a  compiled  language  using  the  integration  routine 
LSODA  [26,41]  that  can  efficiently  accommodate  the  banded 
structure  of  the  system.  When  the  reformer  model  is  simulated 
in  a  native  code  environment  (C  linked  to  LSODA  in  Fortran), 
the  same  simulation  that  took  200  s  using  the  Matlab®  stiff 
system  solver  ode23s  completes  in  about  0.1  s.  However,  ad¬ 
justing  model  parameters,  changing  simulation  protocols,  and 
manipulating  the  simulation  output  are  all  more  cumbersome 
in  a  purely  native  code  environment.  Therefore,  the  external 
C-language  API  for  Matlab®  was  utilized  to  create  an  inter¬ 
face  between  Matlab® /Simulink®  and  the  compiled  reformer 
model/LSODA  solver.  This  allows  for  easier  manipulation  of 
parameters  and  creation  of  simulation  protocols.  Additionally, 
the  simulation  output  can  be  directly  accessed  and  plotted  within 
Matlab®.  The  simulation  of  the  model  through  this  interface 
does  add  some  computational  overhead,  extending  simulation 
time  to  0.7  s,  but  the  utility  and  ease-of-use  outweighs  the  per¬ 
formance  drop  over  native  code  while  still  providing  substantial 
improvement  over  a  pure  Matlab® /Simulink®  implementation. 
With  its  relatively  high  fidelity  and  an  execution  time  of  <  1  s, 
this  model  is  an  excellent  candidate  to  use  for  in  an  optimal 
design  framework. 
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2.1.5.  Model  validation 

The  report  by  Dams  et  al.  [17]  contains  data  from  sev¬ 
eral  experimental  trials,  which  is  useful  for  calibrating  the  re¬ 
former  model  developed  above.  Experiments  were  performed 
in  which  the  external  heating  was  adjusted  alone  or  simultane¬ 
ously  with  feed  flow  adjustments.  Their  test  unit  was  a  stain¬ 
less  steel  reactor  bed  45  cm  long  and  5.2  cm  in  diameter,  cor¬ 
responding  to  an  internal  volume  of  approximately  1  L.  It  was 
filled  with  1.5  kg  of  ICI 3  3-5  A  Cu/Zn0/Al203  low- temperature 
shift  catalyst.  The  outer  wall  was  wrapped  with  a  600  W  elec¬ 
trical  heating  tape  and  covered  with  insulating  material  to  re¬ 
duce  heat  loss  to  the  surroundings.  The  feed  to  the  reformer 
consisted  of  a  3:2  mixture  (volume  basis)  of  methanol  and  wa¬ 
ter,  a  methanol  weight  fraction  of  0.542.  The  liquid  feed  was 
vaporized  and  superheated  prior  to  injection  in  order  to  pre¬ 
vent  condensation  upon  entry  to  the  catalyst  bed.  Since  PEM- 
FCs  operating  on  reformate  are  expected  to  run  at  elevated  an¬ 
ode  pressures,  the  reformer  was  maintained  at  30psig,  even 
though  elevated  pressure  favors  the  production  of  CO  instead 
of  H2. 

Several  steps  were  necessary  to  fit  the  data.  The  experimen¬ 
tal  system  varied  the  power  level  in  the  heating  tape,  whereas 
the  model  assumes  that  heat  is  applied  through  the  variation  of 
wall  temperature.  Since  a  change  in  heater  power  does  not  pro¬ 
duce  an  instantaneous  change  in  wall  temperature,  the  model 
approximates  the  experimental  conditions  by  passing  any  wall 
temperature  changes  through  a  first-order  lag,  with  a  filter  time 
constant  adjusted  to  fit  the  data.  Ultimately,  the  fitting  procedure 
was  an  iterative  process.  First  the  wall  temperature  changes  and 
heat  transfer  coefficient  were  adjusted  to  match  the  temperature 
profile.  Then  the  activation  energy  and  pre-exponential  factor  of 
the  reforming  reaction  were  adjusted  to  align  the  flow  rate  more 
closely.  Next,  the  activation  energy  and  pre-exponential  factor 
of  the  methanol  decomposition  reaction  were  adjusted  to  fit  the 
mole  fraction  of  CO  better.  Altering  the  rates  of  reaction  causes 
the  conversion  to  change,  changing  the  energy  demands.  Thus, 
one  sequentially  adjusts  the  temperature  steps  and  then  the  rate 
parameters,  and  then  iterates  until  convergence  is  achieved.  The 
results  of  that  process  are  shown  in  Fig.  1  with  the  simulation 
represented  by  the  continuous  dashed  line  overlayed  on  the  data 
extracted  from  [17]. 

2.2.  PEMFC  stack 

In  addition  to  the  methanol  reformer,  a  model  is  needed  to 
simulate  a  fuel  cell  stack.  Significant  work  has  been  done  re¬ 
cently  to  create  very  detailed  models  of  individual  fuel  cells 
and  fuel  cell  stacks  that  take  into  account  channel  and  flow-field 
geometries  and  model  micro-  and  meso-scale  transport  phenom¬ 
ena,  for  instance  Shimpalee  et  al.  [45],  Ferng  et  al.  [20],  and  Gu- 
velioglu  and  Stenger  [22] .  Computational  fluid  dynamics  (CFD) 
based  models  are  very  complex,  with  most  requiring  parallel 
methods  to  find  steady- state  solutions,  let  alone  real-time  dy¬ 
namic  simulations.  These  types  of  models  provide  a  great  deal 
of  insight  into  transport-dominated  phenomena  such  as  current 
and  water  distributions  inside  a  single  PEMFC  or  stack;  they  are 
quite  valuable  for  improving  cell  and  stack  flow  channel  designs. 
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Fig.  1 .  Reformer  model  reproducing  the  temperature,  flow  rate,  and  mole  percent 
CO  profile  data  from  Dams  et  al.  [17]. 


However,  CFD  approaches  are  more  detailed  than  what  is  re¬ 
quired  for  balance-of-plant  design  schemes.  These  procedures 
mainly  require  a  correlation  between  load  and  cell/stack  volt¬ 
age,  mass  flows,  and  average  stack  and  stream  temperatures. 
Nevertheless,  these  correlations  are  generally  based  on  preced¬ 
ing  experimental  and  mechanistic  modeling  work. 

At  the  core  of  a  fuel  cell  system  is  the  membrane  electrode 
assembly  (MEA).  Understanding  the  physical  and  chemical  be¬ 
havior  of  the  MEA  is  crucial  to  the  successful  development  of 
a  fuel  cell  model,  but  it  has  not  been  a  simple  task.  At  the  min¬ 
imum,  one  must  consider  multi-component  diffusion  through 
porous  media,  electrochemical  kinetics  such  as  the  reduction  of 
oxygen  in  the  cathode  catalyst  layer,  as  well  as  proton  and  wa¬ 
ter  transport  through  the  electrolyte  membrane.  Consequently,  a 
great  deal  of  effort  has  been  made  to  better  understand  these  phe¬ 
nomena  within  the  last  10-15  years  [9, 10,36,46,25, 19,53,54,43] . 

As  the  technology  and  expertise  regarding  fuel  cells  in¬ 
creased,  models  began  to  appear  that  were  specific  to  particular 
hardware.  Amphlett  et  al.  developed  a  rather  thorough  model 
of  the  Ballard  Mark  IV  PEMFC  using  both  mechanistic  [4]  and 
empirical  [5]  approaches.  The  same  group  was  also  one  of  the 
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first  to  attempt  a  dynamic  PEMFC  model  [8].  Their  efforts  in 
dynamic  modeling  were  later  improved  upon  by  Wohr  et  al. 

[50] .  Using  experience  gained  from  modeling  a  number  of  dif¬ 
ferent  fuel  cells,  Mann  et  al.  [31]  developed  a  generalized  steady 
state  electrochemical  model  of  a  PEMFC  that  can  be  adjusted 
for  individual  cell  parameters  such  as  active  membrane  area  and 
membrane  thickness.  Very  recently,  yet  another  dynamic  fuel 
cell  model  has  been  published  by  Shan  and  Choe  [44] . 

Finally,  the  knowledge  base  developed  for  modeling  individ¬ 
ual  fuel  cells  has  been  combined  to  model  entire  fuel  cell  stacks. 
Much  of  the  behavior  of  the  overall  stack  is  still  related  to  the 
electrochemistry  and  transport  within  the  MEA,  but  greater  dif¬ 
ficulties  are  encountered  in  the  areas  of  heat  and  water  man¬ 
agement.  Since  stack  models  are  typically  used  for  system-level 
applications,  much  of  the  spatial  detail  is  not  needed.  Addition¬ 
ally,  the  most  detailed  MEA  and  individual  fuel  cell  models 
are  often  very  computationally  intensive,  which  would  be  fur¬ 
ther  compounded  as  multiple  cell  models  are  combined  into  a 
stack.  Therefore,  a  number  of  system-level  models  have  been  de¬ 
veloped  in  which  area-  or  volume-averaged  quantities  are  used 
in  addition  to  empirical  correlations  in  order  to  make  the  stack 
models  more  tractable.  An  early  example  of  this  formulation  was 
published  by  Amphlett  et  al.  [6] .  This  approach  was  used  by  Yer- 
ramalla  et  al.  [52]  in  their  dynamic  fuel  cell  model.  Ceraolo  et  al. 
[14]  have  also  developed  a  dynamic  fuel  cell  model  based  on  the 
underlying  physical  phenomena,  but  simplified  enough  to  allow 
integration  into  larger  system-level  simulations.  Recently,  a  dy¬ 
namic  PEMFC  model  was  developed  by  Blackwelder  [11]  that 
can  be  integrated  into  the  University  of  South  Carolina’s  virtual 
test  bed  (VTB)  design  suite.1  Finally,  Xue  et  al.  and  Pathapati 
et  al.  [51,38]  have  developed  a  relatively  simple  system-level 
dynamic  PEMFC  stack  model.  While  this  model  lacks  the  de¬ 
tailed  spatial  information  and  mechanistic  transport  behaviors, 
it  performs  suitably  well  for  overall  material  and  energy  flows 
as  well  as  stack  voltage  with  relatively  minimal  computational 
demand. 

2.2.7.  Model  formulation 

Ultimately,  the  model  presented  simultaneously  by  Xue  et  al. 

[51]  and  Pathapati  et  al.  [38]  was  selected  for  implementation, 
since  it  possessed  the  desired  characteristics:  dynamic  treatment 
of  cell/stack  voltage,  thermal  behavior,  and  material  flows  in  a 
relatively  simple  formulation.  There  are  a  total  of  seven  states 
in  the  model. 

•  Masses  of  N2  and  O2  (idealized  air)  in  the  cathode  as  well  as 

H2  in  the  anode  (3). 

•  Enthalpies/temperatures  of  the  anode  and  cathode  volumes 

and  the  stack  body  (3). 

•  Dynamic  activation  overpotential  Vact  due  to  the  charge 

double-layer  capacitance  effect  (1). 

The  mass  balances  are  very  simple,  since  the  anode  and  cath¬ 
ode  are  considered  to  be  two  well-mixed  volumes  with  a  single 


1  http://www.vtb.engr.sc.edu. 


composition  for  each.  The  balance  equations  are  of  the  general 
form: 


dm/ 

dt 


IP  i,  out 


^Udiff 


(26) 


where  m/jn  and  m/,0ut  are  the  inlet  and  outlet  mass  flow  rates 
of  species  v  through  flow  channels  and  m/,diff  is  the  flow  rate  of 
reactive  species  v  (H2  and  O2)  into  the  MEA.  The  assumption  is 
made  that  the  diffused  species  are  consumed  quickly  once  they 
enter  the  MEA,  so  that  the  rates  of  diffusion  can  be  determined 
stoichiometrically  from  the  current  being  drawn  on  the  stack, 
which  is  given  by  the  following  equation: 


dii,  diff  =  Mi — —N  (27) 

ViT 


Note  that  M\  is  the  molecular  weight  of  species  i,  I  the  cur¬ 
rent  drawn  on  the  stack,  v/  the  stoichiometric  factor  between 
the  number  of  moles  of  species  i  and  the  number  of  moles  of 
electrons  consumed  in  the  reaction  (2  for  H2  and  4  for  O2),  T  is 
Faraday’s  constant,  and  N  is  the  number  of  cells  in  the  stack. 

The  inlet  and  outlet  flow  rates  are  given  by  a  combination  of 
a  simplified  linear  nozzle  correlation  and  the  ideal  gas  law.  The 
nozzle  correlation  assumes  that  the  volumetric  flow  rate  across 
a  valve  is  proportional  to  the  pressure  drop  across  the  valve  as 
the  following  equation: 


V  =  kvAP 


(28) 


Using  the  ideal  gas  law,  the  volumetric  flow  rate  is  converted 
to  a  molar  flow  rate,  which  yields  a  mass  flow  rate  by  use  of  the 
molecular  weight: 

dli, in/out  =  k v  A  P  Mi  (29) 

For  consistency,  the  pressure  P  and  temperature  T  are  always 
assumed  to  be  that  of  the  upstream  side  of  the  valve.  This  flow 
formulation  is  one  of  the  few  differences  from  the  published 
form  of  the  model.  In  [51,38],  the  mass  flow  rate  instead  of 
the  volumetric  flow  rate  was  specified  to  be  proportional  to  the 
pressure  drop.  While  that  may  be  shown  for  pure  gasses  or  fixed- 
composition  gas  mixtures,  this  assumption  seems  rather  dubious. 
Under  ideal  gas  conditions,  the  same  number  of  moles  (equiv¬ 
alently  volume  at  fixed  T ,  P )  should  pass  through  a  valve  for 
a  given  pressure  drop,  relatively  independently  of  the  compo¬ 
sition.  Hence  the  mass  flow  rate  may  vary  considerably  with 
composition,  but  the  volumetric  flow  rate  should  not.  Given  that 
the  composition  of  the  cathode  changes  with  the  electrical  load, 
it  seems  more  appropriate  to  use  a  volume-based  flow  correla¬ 
tion.  The  anode  mass  balance  has  been  modified  similarly  to 
include  CO2  as  an  inert  gas  that  is  part  of  an  idealized  reformate 
mixture. 

The  energy  balances  are  nearly  as  rudimentary.  In  the  anode, 
the  balance  is  written: 

— ^  =  (h A)an?s(Ts  —  Tan)  +  77H2,in 

-  di r2 , out 77h2 , out  -  ^H2, diff  Hr2, diff  (30) 

where  ( hA)m ?s  is  a  lumped  term  containing  the  convective  heat 
transfer  coefficient  and  convective  transfer  area  between  the  an- 
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ode  and  the  stack  body,  and  the  Hn2  terms  are  the  specific  en¬ 
thalpies  of  the  hydrogen  gas  streams  flowing  into  the  anode, 
out  of  the  anode,  and  diffusing  into  the  MEA.  With  the  ideal  gas 
assumption  and  a  constant  heat  capacity  approximation,  the  spe¬ 
cific  enthalpy  for  hydrogen  gas  relative  to  an  arbitrary  reference 
state  is: 

Hu2(T)  =  //H2,ref  +  Cp,H2  (T  ~  Tref)  (31) 

where  Cp, Ho  is  the  constant-pressure  specific  heat  capacity  of 
hydrogen  gas.  A  convenient  choice  for  the  reference  state  is  at 
298.15  K  where  AHfu2  is  0.0  J  kg-1 .  Thus,  the  specific  enthalpy 
simplifies  to: 

HR2(T)  =  CvM2  (T  —  298.15)  (32) 

Note  that  the  enthalpy  does  not  depend  on  pressure,  due  to 
the  ideal  gas  assumption.  The  energy  balance  for  the  stack  body 
has  a  few  additional  terms,  namely  convective  heat  loss  to  the 
surroundings: 

£s,loss  —  (^^)s,room(^room  ^s)?  (33) 


an  equivalent  circuit,  a  differential  equation  can  be  written  for 
the  activation  overpotential,  as  shown  in  the  following  equation: 

dVact  =  /  _  Vact  =  /  /  _  Vkct\ 

d t  C  R,ci  C  C  V  ^act  J  { 

Additionally ,  there  are  ohmic  losses  in  the  stack  predomi¬ 
nantly  related  to  transport  processes.  These  losses  are  found  to 
depend  mainly  on  temperature  and  current  density  and  have  also 
been  previously  determined  to  obey  an  empirical  relationship 
[4,5]  of  the  form 

tfohm  =  ?5+^  +  ?7Z.  (39) 

Finally,  the  cell  voltage  Vceii  is  calculated: 

^cell  —  3^Nernst  Vact  iRohm  (40) 

where  VNemst  is  a  Nemstian  relationship  for  the  equilibrium  cell 
potential  dependent  on  temperature  and  the  partial  pressures 
of  hydrogen  and  oxygen.  The  reader  is  referred  to  [51,38]  and 
references  therein  for  a  more  complete  treatment  of  the  PEMFC 
equations. 


the  heat  from  the  electrochemical  reaction  to  form  water: 

7^s,rxn  —  ^H2,diff  A//rxn  (34) 

(A  Hrxn  normalized  per  kilogram  H2  consumed),  and  energy  loss 
equal  to  the  rate  of  electrical  energy  production: 

^s,elec  —  N  Vcellf  (35) 

where  again,  N  is  the  number  of  cells  in  the  stack,  Vceu  the 
voltage  across  a  single  cell,  and  I  is  the  current  drawn  on  the 
stack.  Finally,  a  generalized  stack  cooling  term  has  been  added 
to  the  original  model  to  allow  for  the  removal  (or  addition)  of 
heat  to  maintain  the  stack  body  at  a  desired  temperature.  It  is  of 
the  form: 

^s,cool  —  Gcool  (36) 

where  (9  cool  becomes  an  additional  adjustable  input  to  the 
model.  This  allows  for  the  construction  of  a  control  loop  to  reg¬ 
ulate  stack  temperature,  although  it  does  not  specify  the  form  or 
geometry  of  any  cooling  solution. 

The  remaining  dynamic  state  of  this  fuel  cell  model  is  the 
activation  overpotential  Vact.  It  is  primarily  due  to  the  slow  ki¬ 
netics  of  oxygen  reduction  in  the  cathode  catalyst  layer.  By  theo¬ 
retically  analyzing  the  physical  phenomena  responsible  for  this 
effect,  the  overpotential  was  determined  to  be  largely  a  function 
of  temperature,  the  current  density  drawn  from  the  cell,  and  the 
effective  concentration  of  oxygen  in  the  cathode  catalyst  layer, 
Cq9.  Using  this  knowledge,  a  steady  state  empirical  correlation 
was  constructed  for  the  activation  overpotential,  given  by  the 
following  equation: 

Oaci  =  §1  +  §2  T  +  §3  T  In  i  +  §4 T  In  c*Ql  (37) 

While  this  correlation  works  well  under  steady  state  con¬ 
ditions,  there  is  an  additional  dynamic  effect  due  to  the  charge 
double  layer  capacitance,  which  arises  due  the  buildup  of  charge 
at  the  membrane/electrode  interface.  By  visualizing  the  stack  as 


2.2.2.  Simulink®  implementation 

Given  its  simplicity,  the  decision  was  made  to  implement  the 
PEMFC  model  directly  within  Matlab®  for  its  ease  of  use  instead 
of  native  code  algorithms  such  as  FSODA  [26,41].  To  provide 
some  performance  enhancement,  the  ODEs  of  the  model  were 
implemented  as  a  compiled  C-MEX  S -function.  A  variety  of  pro¬ 
cess  inputs  and  model  parameters  were  left  exposed  externally 
in  order  to  be  adjusted  rapidly  within  the  Simulink®  environ¬ 
ment  without  recompiling.  In  addition  to  the  seven  states  in  the 
model,  a  number  of  other  useful  quantities  were  output  from  the 
model,  including  anode/cathode  pressures,  stack  voltage,  and 
inlet/outlet  mass  flow  rates. 

The  model  requires  volumes  for  the  cathode  and  anode  chan¬ 
nels  in  order  to  compute  the  respective  mass  balances.  However, 
for  the  purpose  of  performance  and  power  scaling,  the  mem¬ 
brane  area  is  the  critical  parameter.  Therefore,  considering  the 
envisioned  use  for  an  optimal  design  problem,  the  worksheet 
specifies  a  cell  membrane  area  and  number  of  cells  in  the  stack 
and  then  assumes  constant  volume/total  area  ratios  to  calculate 
the  volumes  of  the  anode  and  cathode.  For  the  test  cases  pre¬ 
sented  in  the  work  of  Xue  et  al.  and  Pathapati  et  al.  [51,38],  the 
cathode  was  specified  with  a  volume  twice  that  of  the  anode. 
That  ratio  is  preserved  in  this  worksheet.  Based  on  rough  ge¬ 
ometrical  arguments,  the  ratio  for  the  anode  was  calculated  to 
be  0.2  and  0.4  cm3  cm-2  for  the  cathode.  Similarly,  the  weight 
of  the  stack  is  computed  based  on  the  cell  membrane  area  and 
number  of  cells.  Additionally,  the  current  density  drawn  from 
a  stack  is  more  informative  than  the  actual  current,  so  the  load 
on  the  stack  is  specified  in  terms  of  current  density  and  then 
multiplied  by  the  membrane  area  to  determine  the  actual  loads 
on  the  PEMFC  stack. 

2.3.  Heat  exchanger 

Due  to  the  differing  operating  temperatures  between  different 
components  within  the  system,  it  will  be  necessary  to  heat  and 
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cool  a  number  of  different  process  streams.  In  order  to  improve 
efficiency,  it  is  possible  to  use  excess  heat  from  a  hot  stream  such 
as  the  reformate  to  heat  other  streams  such  as  the  cathode  air  feed 
to  the  reformer  stack.  Cross-flow  heat  exchanger  geometries  are 
common  for  gas-gas  heat  exchange,  so  a  simple  steady  state 
heat  exchanger  model  is  used. 

Assuming  uniform  velocity  profiles  and  constant  physical 
properties,  the  governing  equations  for  the  heat  exchanger  can 
be  made  dimensionless  with  the  solution  depending  solely  on  the 
Nusselt  numbers  on  the  hot  and  cold  sides  of  the  heat  exchanger. 
The  dimensionless  problem  was  discretized  using  finite  differ¬ 
ences  and  solved  off-line  for  combinations  of  Nusselt  numbers 
ranging  from  ICC2  to  102.  The  dimensionless  average  hot  and 
cold  outlet  temperature  were  tabulated  for  each  case,  allowing 
the  heat  exchanger  to  be  implemented  in  simulation  as  a  look-up 
table. 

2.4.  Auxiliary  power  source 

While  fuel  cell  stacks  are  typically  able  to  respond  to  large 
variations  in  load  or  power  quickly,  the  time  scales  of  the  fuel 
processing  system  are  generally  at  least  two  orders  of  magnitude 
slower.  Therefore,  the  rate  at  which  the  load  can  be  increased  on 
a  PEMFC  stack  coupled  to  a  fuel  processor  is  limited  by  the  rate 
at  which  hydrogen  production  can  be  increased.  Consequently, 
some  form  of  auxiliary  power  will  be  required  to  ensure  satis¬ 
factory  vehicle  performance  under  normal  operating  conditions. 

The  selection  of  a  specific  auxiliary  power  source,  such  as 
lead-acid  batteries  or  supercapacitors,  is  beyond  the  current 
scope  of  this  work.  Therefore  the  system-level  model  developed 
treats  the  auxiliary  power  source  as  a  black-box  component.  It  is 
assumed  that  this  power  source  operates  at  a  constant  bus  voltage 
-120V  for  this  work  -  and  is  instantaneously  able  to  provide 
enough  current  to  provide  the  difference  between  the  desired 
power  and  the  PEMFC  stack  power.  When  the  PEMFC  is  pro¬ 
ducing  more  power  than  is  required,  the  excess  is  then  used  to 
recharge  the  auxiliary  power  unit.  Furthermore,  there  is  assumed 
to  be  a  5%  loss  in  each  direction  (charging/discharging)  as  part 
of  the  power  conversion  process.  One  is  then  able  to  observe 
the  amount  of  auxiliary  capacity  needed  for  a  particular  system 
configuration  and  power  profile.  The  amount  of  capacity  and  the 
rate  at  which  it  is  drawn  can  then  be  used  to  determine  what  type 
of  auxiliary  power  source  is  most  suitable  for  the  application. 

3.  System-level  integration  and  control 

The  components  of  the  overall  system  model  were  connected 
as  shown  in  Fig.  2  using  Matlab®  and  Simulink®.  This  diagram 
depicts  an  open-loop  configuration.  By  specifying  material  flow 
rates,  operating  temperatures  and  pressures,  and  the  load  on  the 
stack,  outputs  such  as  stack  power  and  hydrogen  production  can 
be  observed.  However,  this  is  not  the  desired  mode  of  operation. 
For  vehicular  operation,  it  is  necessary  to  specify  a  desired  power 
level  and  then  have  controllers  on  the  system  adjust  stack  load 
and  hydrogen  production  to  meet  that  goal.  For  this  work,  an 
idealized  driving  power  profile  was  developed  and  is  shown  in 
Fig.  3.  Although  a  true  power  profile  would  vary  much  more 
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Fig.  2.  Integrated  system  layout  for  the  reformer/PEMFC  stack  system. 
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Fig.  3.  Performance  of  system  following  the  designated  power  profile.  The 
parameter  set  for  this  run  is  {Lref  =  0.8,  Rre f  =  0.07,  Amem  =  900,  TsP)ref  = 
550,  rSp5Stack  =  355}.  Upper:  PEMFC  power  compared  to  the  desired  system 
power.  The  auxiliary  power  at  any  time  is  the  difference  between  the  desired 
power  and  the  PEMFC  power.  Lower:  relative  change  in  capacity  of  the  auxiliary 
power  system  at  the  operating  voltage  of  120  V. 


frequently,  this  profile  is  of  the  correct  magnitude  and  attempts 
to  account  for  warm-up,  acceleration,  braking,  and  cruising. 

Several  control  loops  are  necessary  to  achieve  satisfactory 
closed-loop  operation.  Two  low-level  regulatory  proportional- 
integral  (PI)  feedback  control  loops  provide  temperature  control 
on  the  reformer  and  the  PEMFC  stack.  The  reformate  tempera¬ 
ture  can  be  controlled  by  the  manipulation  of  the  reformer  wall 
temperature.  In  actual  practice,  the  wall  temperature  is  not  the 
true  manipulated  variable;  it  would  likely  be  a  valve  controlling 
fuel  flow  to  a  burner  element.  Similarly,  the  temperature  of  the 
PEMFC  is  maintained  by  adjusting  the  fictitious  (2  cool  term.  In 
practice,  the  temperature  would  be  regulated  by  adjusting  the 
flow  rate  of  a  cooling  liquid.  Since  the  cooling  fluid  could  never 
be  used  to  heat  the  stack,  the  output  of  the  PI  controller  passes 
through  a  saturation  block  to  disallow  values  of  QCoo\  that  add 
heat  to  the  stack.  Tunings  for  the  controller  gains  Kc  and  the 
integral  time  constant  x\  are  given  in  Table  2. 

Two  additional  loops  have  the  greatest  impact  on  the  perfor¬ 
mance  of  the  closed-loop  system.  The  first  is  the  regulation  of  the 
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methanol  flow  rate  into  the  reformer  system.  The  second  is  the 
regulation  of  the  load  on  the  PEMFC.  The  manipulation  of  both 
quantities  is  based  on  the  difference  between  the  desired  system 
power  and  the  PEMFC  stack  power.  For  the  sake  of  simplic¬ 
ity,  these  controllers  were  also  modified  forms  of  PI  controllers. 
Ultimately,  the  implementation  was  more  nuanced  for  several 
reasons.  First,  the  power  response  of  the  PEMFC  displays  direct 
feedthrough  since, 


^stack  —  i  Amem  Vstack- 


(41) 


This  direct  feedthrough  seemed  to  induce  instabilities  ex¬ 
cept  at  extremely  low  controller  gains  when  using  a  standard 
PI  formulation,  most  likely  due  to  the  large  gain  on  the  di¬ 
rect  feedthrough.  In  order  to  mitigate  the  impact  of  the  direct 
feedthrough  (and  equivalently  step  changes  in  the  setpoint),  the 
error  signal  (Psp  —  Pstack)  is  first  passed  through  a  first  order 
filter  before  reaching  the  controller.  Thus,  the  control  law  is  of 
the  form: 

m  =  Kc  ( 1  +  —)  - — (PSp  (s)  -  PstackW),  (42) 

V  T\S  J  TfS+1 

where  s  is  the  Faplace  variable.  The  control  structure  is  identical 
for  the  methanol  flow  controller.  For  the  power  controller,  a  filter 
time  constant  of  0. 1  min  was  selected,  while  a  larger  value  of 
2  min  was  used  for  the  methanol  flow  rate,  given  the  slower 
response  of  that  component.  Additionally,  the  power  setpoint 
Psp  was  more  complicated  than  simply  the  desired  power  output 
at  a  given  time.  It  was  desired  to  operate  the  system  so  that  the 
auxiliary  power  source  remained  as  close  to  a  reference  charge 
level  as  possible.  By  integrating  the  load  on  the  auxiliary  power 
source,  it  was  possible  to  track  the  relative  change  in  capacity  of 
that  system.  This  auxiliary  capacity  level  was  then  used  to  adjust 
the  power  setpoint  according  to  the  following  linear  feedback 
law: 


ftp  —  ^desired  T  ^"aux  Caux.  (43) 

For  power  measured  in  Watts  and  capacity  in  A  h,  a  feedback 
proportionality  constant  Kaux  of  — 1000  W  (A  h)-1  provided  sat¬ 
isfactory  response  without  substantial  oscillation  in  the  auxiliary 
capacity  level.  Finally,  actuator  limits  were  placed  on  the  con¬ 
trol  actions  in  accordance  with  physical  limits.  The  methanol 
flow  rate  and  load  on  the  fuel  cell  stack  cannot  be  lower  than 
0  F  min- 1  and  0  A  cm-2 .  The  maximum  methanol  flow  rate  was 
established  as  1  F  min-1 ,  which  was  somewhat  arbitrary,  but  in¬ 
tentionally  left  rather  high.  Finally,  the  MEAs  in  the  PEMFC 
stack  can  be  damaged  if  the  unit  is  starved  of  hydrogen,  so  there 
is  a  dynamic  upper  limit  on  PEMFC  load,  preventing  the  load 
from  exceeding  that  needed  to  consume  85%  of  the  hydrogen 


Fig.  4.  Control  actions  taken  by  the  major  system  controllers  corresponding  to 
the  desired  power  trajectory  shown  in  Fig.  3.  The  parameter  set  for  this  run 
is  {Lref  =  0.8,  Rrcf  —  0.07,  Amem  =  900,  TSpjref  =  550,  7spjStack  =  355}.  Up¬ 
per:  liquid  methanol  flow  rate  into  the  reformer.  The  liquid  water  flow  rate  into 
the  reformer  is  fixed  at  a  2:3  ratio  to  the  methanol  flow  rate.  Lower:  load  current 
density  (solid)  on  the  PEMFC  stack  as  a  function  of  time.  The  dynamic  upper 
constraint  (dashed)  is  the  current  density  which  would  use  85%  of  the  hydrogen 
output  from  the  reformer. 

being  produced  by  the  reformer.  Representative  control  action 
for  these  two  controllers  is  shown  in  Fig.  4. 

4.  Results  and  discussion 

Simulation  studies  were  performed  to  gauge  the  effects  of 
various  design  and  operational  parameters  on  the  overall  perfor¬ 
mance  of  the  system.  There  are  a  number  of  parameters  that  may 
influence  operation,  only  five  were  chosen  for  consideration  in 
this  study.  Three  were  design  parameter  related  to  the  physical 
structure  of  the  system:  reformer  length,  reformer  radius,  and  ac¬ 
tive  membrane  area  for  each  cell  of  the  PEMFC  stack.  The  two 
additional  parameters  were  operational  in  nature:  the  reformer 
outlet  temperature  setpoint  and  the  PEMFC  stack  body  temper¬ 
ature  setpoint.  Four  values  of  each  parameter  were  considered 
and  are  listed  in  Table  3.  Values  for  the  other  parameters  used 
in  the  system  model  are  contained  in  Table  4. 

Given  5  parameters,  each  with  4  possible  values  allows  for 
1024  different  system  configurations  to  be  studied.  For  each  pa- 


Table  2 


Control  loop  tuning  parameters 


Loop 

Kc 

Ti  (min) 

Reformer  temperature 

IKK-1 

10 

PEMFC  temperature 

5kJK-1 

10 

PEMFC  load 

1(T4  (A  cm-2)  W-1 

1 

MeOH  flow 

10-5  (L  min-1)  W-1 

2 

Table  3 


Parameter  values  studied  in  the  sensitivity  analysis 


Parameter 

1 

2 

3 

4 

Reformer  length  (m) 

0.600 

0.80 

1.00 

1.20 

Reformer  radius  (m) 

0.025 

0.04 

0.07 

0.10 

PEMFC  membrane  area  (cm2) 

800 

900 

1000 

1200 

Reformer  temperature  setpoint  (K) 

500 

550 

575 

600 

PEMFC  temperature  setpoint  (K) 

335 

345 

355 

365 
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rameter  combination,  the  system  is  simulated  trying  to  reproduce 
the  hypothetical  power  profile  shown  in  Fig.  3.  No  effort  is  made 
to  tune  the  control  loops  for  each  configuration.  Consequently, 
the  system  response  for  each  configuration  may  not  be  ideal.  In 
fact,  results  were  not  obtained  for  several  configurations,  due  to 
numerical  difficulties,  possibly  resulting  from  controller  tun¬ 
ings  being  unstable  at  the  given  configuration.  Nevertheless, 
system  performance  was  analyzed  based  on  factors  including 
total  methanol  usage,  hydrogen  utilization  percentage,  auxiliary 
power  capacity  range,  and  a  gross  efficiency  measure  for  all 
configurations  for  which  results  are  available. 

Fig.  5  is  a  histogram  of  the  net  methanol  usage  for  the  vari¬ 
ous  system  parameter  combinations  over  the  one  hour  reference 
power  trajectory.  Note  that  the  distribution  appears  bimodal  with 
a  large  normal-looking  distribution  centered  around  1 1  L  and  a 
smaller  peak  distributed  between  about  6.5  and  8.0  L.  By  ana¬ 
lyzing  the  parameter  combinations  that  were  found  in  the  sec¬ 
ond  peak,  it  was  found  that  all  were  at  the  lowest  value  of  the 
reformer  setpoint  temperature:  500  K.  Intuition  indicates  that 
lower  methanol  usage  should  be  a  positive  characteristic  for  a 
design,  but  system  design  must  consider  a  variety  of  factors. 
Before  reaching  any  conclusions,  it  is  necessary  to  examine  the 
performance  of  these  parameter  configurations  in  the  other  per¬ 
formance  categories.  The  next  category  is  hydrogen  utilization. 
For  this  analysis,  hydrogen  utilization  is  reported  as  the  per¬ 
centage  of  hydrogen  produced  by  the  reformer  that  is  consumed 
by  a  single  pass  through  the  PEMFC  stack.  In  practice,  it  may 
be  possible  to  run  anode  exhaust  into  a  recycle  loop  or  at  least 
flare  the  remaining  hydrogen  to  capture  its  energy  as  heat,  but 
these  options  are  not  considered  here.  The  distribution  of  hy¬ 
drogen  utilization  over  the  parameter  combinations  is  shown  in 
Fig.  6.  Utilization  is  distributed  widely  between  about  55%  and 
78%  with  a  large  grouping  between  55%  and  60%.  Furthermore, 
it  was  found  that  the  parameter  combinations  that  were  in  the 
smaller  peak  for  methanol  utilization  were  also  the  combinations 
that  had  greater  than  70%  hydrogen  utilization.  Once  again,  it 
would  appear  that  high  hydrogen  utilization  would  be  a  positive 
design  characteristic. 

However,  the  picture  changes  dramatically  when  auxiliary 
power  capacity  is  considered.  Fig.  7  illustrates  the  distribution 
of  auxiliary  power  capacity  for  the  various  combinations.  Once 

Table  4 


Other  system  operational  parameters 


Parameter 

Value 

Number  of  cells  (#) 

120 

Aux.  power  voltage  (V) 

120 

PEMFC  cathode  pressure  (atm) 

3.0 

Reformer  pressure  (atm) 

2.0 

Reformer  heat  transfer  coefficient  (J  m-2  s  K) 

50.0 

Liquid  MeOH/H20  ratio  (v/v) 

3  :  2 

HX  height  (m) 

0.10 

HX  cold  length  (m) 

0.15 

HX  hot  length  (m) 

0.10 

HX  channels/side  (#) 

5 

HX  hot  channel  Cp  (J  mol- 1  K) 

36.0 

HX  cold  channel  Cp  (J  mol-1  K) 

28.0 

HX  heat  transfer  coefficient  (J  m-2  s  K) 

50.0 

120 


Net  Methanol  Usage  (L) 


Fig.  5.  Histogram  of  methanol  usage  for  the  various  parameter  configurations 
following  the  reference  power  trajectory  shown  in  Fig.  3.  Net  usage  is  computed 
assuming  that  95%  of  unreacted  methanol  leaving  the  reformer  can  be  recaptured 
and  reused. 

again,  there  is  a  very  wide  range  of  results  between  about  10 
and  120  Ah  at  120  V  operating  voltage  with  a  significant  num¬ 
ber  (58%)  below  20  A  h.  The  amount  of  surge  power  needed  in 
the  system  is  very  critical  to  the  overall  design,  since  it  strongly 
influences  both  the  weight  and  the  cost  of  the  system.  This  bodes 
poorly  for  the  low  methanol  consumption/high  hydrogen  utiliza¬ 
tion  parameter  combinations,  since  they  were  the  ones  which 
required  80  A  h  and  higher  of  auxiliary  capacity.  This  represents 
at  least  a  four-  to  six-fold  increase  over  the  majority  of  the  other 
configurations,  and  suggests  that  these  particular  designs  are  in¬ 
adequate. 

Finally,  the  histogram  of  gross  efficiencies  is  presented  in  Fig. 
8.  For  this  study,  the  gross  efficiency  was  considered  to  be  the 
sum  of  the  electrical  output  of  the  PEMFC  stack  and  the  change 
in  the  electrical  energy  of  the  auxiliary  power  divided  by  the 
higher  heating  value  of  the  methanol  consumed  (17,990  kJ  L-1). 

300 1 - 1 - 1 - 1 - 1 - 1 - 


Net  H2  Utilization  (%) 


Fig.  6.  Histogram  of  hydrogen  utilization  for  the  various  parameter  configura¬ 
tions  following  the  reference  power  trajectory  shown  in  Fig.  3. 
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Aux.  Power  Capacity  (A-hr) 


Fig.  7.  Histogram  of  the  necessary  sizing  of  the  auxiliary  power  capacity  needed 
for  the  various  parameter  configurations  to  follow  the  reference  power  trajectory 
shown  in  Fig.  3. 

There  are  a  number  of  factors  not  considered  in  this  measure, 
including  the  heat  of  vaporization  for  the  methanol  and  water 
feed  to  the  reformer,  the  external  heat  to  drive  the  reactions  in 
the  reformer,  the  heating  value  of  the  hydrogen  in  the  PEMFC 
anode  exhaust,  and  other  ancillary  costs  such  as  pumping  costs 
and  work  to  cool  the  PEMFC  stack.  Consequently,  these  num¬ 
bers  should  not  be  viewed  as  an  absolute  quantitative  measure, 
but  rather  a  relative  scale  to  compare  performance  among  var¬ 
ious  parameter  combinations.  Most  of  the  results  are  normally 
distributed  around  an  efficiency  value  of  about  36%,  but  again 
there  is  an  extended  lower  tail.  As  expected,  the  low  methanol 
use  combinations  comprise  the  lower  tail  in  the  range  of  20-30%. 

More  insight  is  gained  by  examining  the  dynamic  profile  of 
these  low  methanol  consumption  cases.  It  is  apparent  that  at  low 
reformer  operating  temperatures,  the  conversion  of  methanol  to 
hydrogen  is  very  low,  regardless  of  the  methanol  flow  rate,  result¬ 
ing  in  the  low  net  methanol  usage.  Thus,  due  to  the  control  con- 


140 


Gross  Efficiency  (%) 


Fig.  8.  Histogram  of  the  gross  efficiency  measure  for  the  various  parameter 
configurations  when  following  the  reference  power  trajectory  shown  in  Fig.  3. 


straints  based  on  hydrogen  flow,  the  load  on  the  PEMFC  stack 
is  kept  low  and  cannot  meet  the  power  demand  at  most  times. 
Even  when  using  most  of  the  available  hydrogen,  the  system  is 
forced  to  draw  heavily  on  the  auxiliary  power  unit,  which  is  why 
these  configurations  required  much  more  auxiliary  capacity.  Fi¬ 
nally,  even  at  the  end  of  duty  cycle,  these  configurations  were 
not  able  to  recharge  the  auxiliary  capacity  to  its  initial  level.  This 
net  loss  of  energy  from  the  auxiliary  power  unit  results  in  the 
lowered  efficiencies  of  these  designs,  and  further  indicates  that 
they  would  not  be  good  choices  for  a  vehicular  power  system. 

5.  Future  work 

This  work  begins  to  quantify  the  performance  tradeoffs  that 
result  from  design  and  operational  decisions  for  the  methanol 
reformer/PEMFC  system.  There  are  a  number  of  avenues  to  ex¬ 
plore  this  work  further.  System  performance  can  be  studied  for 
different  equipment  configurations.  Model  improvements  and 
additions  can  be  incorporated.  Currently,  components  such  as  the 
reformer  feed  evaporator  and  any  reformate  processing  such  as 
preferential  oxidation  are  not  considered,  and  while  they  would 
be  expected  to  significantly  impact  the  dynamics  of  the  pro¬ 
cess,  would  give  a  more  accurate  estimate  of  fuel  and  energy 
efficiency  for  the  system.  Additionally,  the  representative  power 
profile  used  to  drive  the  system  is  very  simplistic.  Eventually  it 
would  be  useful  to  consider  a  power  profile  that  is  more  realis¬ 
tic.  By  making  certain  assumptions  about  the  mass  and  drag  of 
the  vehicle,  it  is  possible  to  take  driving  schedules  such  as  the 
federal  urban  driving  schedule  (FUDS)  or  the  federal  highway 
driving  schedule  (FHDS),  which  specify  velocities,  and  calcu¬ 
late  power  profiles.  Since  the  mass  of  the  vehicle  depends  on 
the  size  of  the  main  components,  the  power  profile  should  be 
adjusted  accordingly  for  each  configuration.  Ogden  et  al.  have 
a  clear  demonstration  of  this  approach  [37]. 

Furthermore,  the  sensitivity  studies  were  only  conducted  with 
five  parameters  and  four  values  for  each  parameter.  It  is  therefore 
unlikely  that  any  of  the  exact  parameter  combinations  explored 
here  will  be  optimal  in  any  of  the  performance  metrics.  Given 
that  there  are  a  number  of  different  objectives  to  consider  when 
designing  such  a  power  plant  including  cost,  weight,  size,  and 
efficiency,  this  model  could  be  used  as  part  of  an  multiobjective 
optimization  procedure.  One  approach  would  be  to  develop  a 
single  composite  objective  function,  which  is  a  weighted  com¬ 
bination  of  the  objectives  under  consideration,  similar  to: 

®  =  Gost^cost  T  Gize^size  T  Gnass^mass  Gff<Peff-  (44) 

In  order  to  achieve  this,  one  would  also  need  to  establish  cor¬ 
relations  between  the  easily  calculated  performance  measures 
such  as  auxiliary  power  capacity  and  the  design  objectives  such 
as  monetary  cost,  system  size  (volume),  and  system  mass.  Also, 
there  can  be  considerable  difficulty  in  choosing  appropriate  val¬ 
ues  for  the  weighting  factors  G  that  produce  a  satisfactory  result. 
It  may  be  possible  to  use  techniques  from  Pareto  multiobjective 
optimization  theory  to  obtain  these  values. 

Additionally,  any  optimization-based  approaches  that  include 
dynamic  simulations  as  part  of  the  objective  function  raise  legit¬ 
imate  concerns  about  the  computational  demands  of  the  method. 
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However,  any  optimization  techniques  that  use  finite  differenc¬ 
ing  schemes  to  obtain  gradients  or  otherwise  need  to  evaluate 
multiple  independent  points  in  the  solution  space  can  benefit 
from  parallel  computing  techniques.  For  instance,  for  a  function 
of  N  different  variables,  N  +  1  points  need  to  be  evaluated  to 
compute  the  gradient  at  a  single  point;  once  these  points  are  cho¬ 
sen,  each  can  be  evaluated  independently  without  knowledge  of 
the  values  at  any  other  points.  Thus,  this  sort  of  operation  is  well- 
suited  for  parallelization,  and  can  make  a  significant  impact  on 
the  time  necessary  for  optimization  calculations  [13,24]. 

Clearly,  the  development  of  liquid-fueled  vehicular  fuel  cell 
power  plants  is  a  very  active  area  for  ongoing  research,  both 
in  the  fundamental  theoretical  understanding  of  the  processes 
involved  and  the  application  of  techniques  for  design  and  inte¬ 
gration. 

6.  Conclusions 

This  work  describes  the  development  and  implementation 
of  a  system  level  model  of  a  vehicular  reformer/PEMFC  stack 
power  system.  Particular  emphasis  has  been  placed  on  the  in¬ 
terconnections  between  sub-units  and  the  formulation  of  con¬ 
trol  loops  to  operate  the  system  to  meet  desired  power  targets. 
With  a  particular  hardware  configuration,  a  sensitivity  study  was 
conducted  to  determine  the  influence  of  various  design  and  op¬ 
erating  parameters  on  system  performance.  Specifically,  it  was 
shown  that  the  operation  of  the  reformer  at  low  temperatures 
gives  particularly  inferior  performance.  Recommendations  for 
a  particular  design  have  not  yet  been  made,  since  they  are  very 
sensitive  to  the  relative  importance  of  the  performance  metrics, 
but  doing  so  remains  a  future  objective.  Nevertheless,  a  useful 
model  framework  has  been  developed  for  the  dynamic  analysis 
of  a  complex  integrated  system. 
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